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A n'ethod, based on the n^aximuni n^edian hourly value of Ji^F^ n?onth-by-n>onth as 
observed at Washington, D.C., is given for predicting a smoothed annual sunspot ruimber 
immediately at the close of a given month, centered on the month in question. Regression 
equations and standard errors are given. It should be a useful supplement to the McNish- 
Lincoln method of prediction of sunspot numbers, particularly during the first two years 
of the rising part of a sunspot cycle. It is capable of use with the observations from any 
ionosphere sounding station operated in a consistent manner over a period of years, prefer- 
ablv during at least two solar cvcles. 



1. Introduction 

One of tho widely acccptcnl methods us(h1 in the 
prediction of sunspot numbers foi" use in radio 
frequency [)r(Hlictions is the A IcNish- Lincoln juethod 
[1949]. This method involves the use of smoothed 
annual sunspot numbers corrected by departures one 
year earlier than the time in question. It predicts 
month by month the ''best" prediction for 1 year 
ahc^ad. 

The use of a 12-month running average suffers 
of course from the fact that to obtain the sunspot 
number a])})licable to a certain month, 6 months in 
the future nuist be observed. Wlien, however, the 
authors extended the technique to tlie prediction of 
3-month mean sunspot numbers, lower overall 
reliability residted. The nu^thod outlined here has 
the advantage^ that, inmiediately upon the comple- 
tion of a month, a prediction of the 12-month running- 
average sunspot number centered on that month can 
be made. Of considerable importance is the fact 
that the rms deviation of prediction would seem to 
be of about the same order of magnitude as those 
given in a second paper [McNish and Lincoln, 1954] 
in which the authors examined the reliability of 
their method of prediction and applied it to cycle 18. 

In the proposed method, the highest monthly 
median value of /o^2 for any hour, as observed at 
Washington, D.C., is used, in a set of 12 regression 
equations, one for each month, to obtain an estimate 
of the 12-month running average sunspot number 
for the month that has just ended. It is believed 
that this choice of maximum /oi^2, rather than /0F2 
for a fixed time of day, is of aid in minimizing 
seasonal effects. 

The idea of using/oF2 as a measure of solar activity 
is not new, but the methods so far proposed have 
been more elaborate, involving for instance such 
things as 5-hoiirly averages centered on noon, and 
averaging of results from several stations [Fliillips, 
1947; Mimus, 1955]. Also nmch work has been 



done and is continuing on other possibh* hidices 
[Kundu and Denisse, 195S; Minnis and Bazzard, 
1959]. The fact that a reasonable answer can be 
obtained to the problem of predicting solar activity 
by the simple method outlined here, gives nnison 
for hope that the nuich more elaborate nuMhods 
presently unch^r study will eventually resuh i!i a Inghly 
satisfactory inckvx obtained from some combina- 
tion of physical (juantities. At the same time, it is 
distinctly wortliwhile to have some way of coimect- 
ing future ionospheric indices with the long time- 
series of sunspot niunbers now available, on wliich 
so many relationships linve previously b(HMi [)ased. 
The method of this pajXM* is intemhHl to show 
what can be (h)ne almost immediately using data 
from but one ionosphere station (Washington) which 
has been operated in a consistent fashion over a 
long period of time. The initial impetus for this 
study was the sudden drop in solar activity in 
January 1961 as represented by the provisional 
Zurich monthly sunspot mimber of 53.5, together 
with the desire to investigate what this sudden 
drop might mean in terms of future predictions of 
critical frequencies. 

2. Correlation Coefficients 

Various results are given in table 1. It will be 
noted that, in general, the coefficients in cohnnn 
A exceed those in column B, and those in cohnnn 
E exceed those in column D. (However, tlu?se 
observed tendencies by themselves are not statisti- 
calty significant.) Column C is filled out only 
for January and May, the months for which th(^ 
highest and lowest coefficients were obtained hi 
column A. The omission of numbers for the maxi- 
mum period 1957 to 1959 (col. C) and the re- 
duction of high numbers in this period (col. E) 
gave correlations which, when taken with thosc^ in 
columns A and D, form an interc^sting scM'ies of 
coefficients for these montlis. 
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Table 1. Correlation coefficients, niaximu7n hourly ynedian 
/0F2 {Washington) versus Zurich sunspot mmiher for various 
periods; A, C, D, versus 12-month running average sunspot 
number centered on the month in question; B, versus monthly 
sunspot number; E, same as D, but with sunspot numbers 
^150 used as 160 



Jan_ 
Feb. 
Mar. 
Apr_ 
May 
June 

July. 
Aug. 
Sep.. 
Oct. 
Nov. 
Dec- 



1945-1960 



97 



1945-1960 



0.97 
.94 
97 
92 
91 



1941-1956, 1960 



0.985 



1941-1960 



.97 
.97 
.93 
.96 

.96 
.95 
.98 
.97 



The use of a value of 150 for all running average 
sunspot numbers over 150, produced an equal or 
better correlation, except for January. This is in 
agreement with previously noted results [Bennington, 
1958; Ostrow, 1959] and with CRPL predictions 
practice during the recent maximum.^ In January, 
at the recent maximum period, increased solar 
activity did seem to have a marked effect in raising 
the criticals at Washington to a very high value 
in the middle of the day, higher than CRPL pre- 
dictions based on sunspot number 150. Also it 

1 CRPL Series D, Basic Radio Propagation Predictions, U.S. Government 
Printing Office, Washington 25, D.C., Nos. 146-173. Predictions for Jan. 1957- 
Apr. 1959, issued Oct. 1956-Jan. 1959. 
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Figure 1. Scatter of points and regression line for March 
{tables 1 and 2, col. E). 



will be noted that correlations are somewhat higher 
for the winter months, a not unexpected result, as it 
has long been noted at Washington that the range 
of variation in /oi^2 over the solar cycle, is much 
ereater for winter months than for summer months. 
While it might be argued that the apparent seasonal 
effect was due to the way the sunspot numbers 
happened to fall, e.g., if there had been four num- 
bers for January around 200 and none for June 
greater than 150, it should be remarked that during 
the recent maximum (1956 to 1959) each montli 
had one 12-month running average sunspot number 
>187, and each month had three >150. Column 
B is included to show that correlations can be very 
liigh between monthly sunspot numbers and critical 
frequency. Figure 1 illustrates tlie relationship 
between / and E leading to the correlation coef- 
ficient 0.97 in column E, table 1, for March. 

3. Regression Equations 

Table 2 gives the linear regression equations for 
each month of the year, based on the coefficients in 
column E of table 1. In addition, equations are 
given for October through February using the 
coefficients of column D, which are greater than or 
equal to the corresponding coefficients in column E. 

Table 2. Regression equations; R=12-rnonth running average 
sunspot number, f= maximum hourly median fQF2{Mcls) ; 
columns D, E correspond to D, E in table 1 



Jan - - 
Feb._ 
Mar.. 
Apr_. 
May. 
June. 

July.. 
Aug.. 
Sep-- 
Oct... 
Nov.- 
Dec. 



P = 83.2+20.5(/-9.6) 
/?=83.4+21.7(/-9.6) 



./?= 84.2+22.9 (/-lO.O) 
/?=84.2+21.4(/-10.3) 
/?=84.0+22.9(/-10.0) 



^ = 78.4+17.8(/-9.6) 
.P = 78.4+19.0(/-9.6) 
/? = 78.4+18.8(/-9.2) 
/? = 78.8+24.1 (/-8.2) 
P = 79.2+39.4(/-7.3) 
/? = 79.5+56.0(/-6.8) 

/? = 79.8+53.4(/-6.4) 
/? =80.0+44.1 (/-6.8) 
P = 79.9+26.3(/-8.1) 
/? = 79.8+20.3 (/-lO.O) 
/? = 79.6+18.8(/-10.3) 
/? = 79.4+20.2(/-10.0) 



4. Standard Error of an Estimated Sunspot 
Number 

Table 3 gives the low and high values of/ (maxi- 
mum hourly median /o^2 in Mc/s) month by month 
over the 20-year period 1941 to 1960, as well as the 
means, which, of course, are repeated from table 

2. Table 4 gives the standard error of an estimated 
sunspot number, R, for each value of/ given in table 

3. These are computed from the formula 



<^/e(i 



-Sr 



■V 



l+i- 

n 



(f~fr 



j:(fi-fy 



[Snedecor, 194(3] 



wdiere n is the number of observed pairs (fi, Ri), 
f=y^,filfi, SR.fis the sample standard deviation about 
the regression line, and / is the predicting value of 
/o/^2 as before. They may be compared with the 
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Table :\. Range of values of f (Mc/s), Washin (/ton, 1941-1960 





Jan. 


Feb. 


Mar. 


Apr. 


May 


June 


July 


Aug. 


Sep. 


Oct. 


Nov. 


Dec. 


Low ._ . _._ _ - --. - 


5.9 
9.0 
15.5 


5.7 
9.6 
14.0 


5.4 
9.2 
13.2 


5.3 

8.2 
11.2 


5.2 

7.3 

10.0 


5.3 

(i. 8 

8.0 


4.8 
6.4 

7.7 


5.0 
6.8 
8.6 


5.2 
8.1 
11.0 


6.2 
10.0 
14.2 


6.5 
10.3 

14.8 


6.2 


Mean . .- 


10.0 


High 


14.0 











Table 4. 


Standard error o 


f an 


estimated R for 


various values 


of f {see table 3) 










Jan. 
Col. D, 
table 2 


Jan. 
Col. E, 
table 2 


Feb. 
Col. D, 
table 2 


Feb. 
Col. E, 
table 2 


Col. E, table 2 


Oct. 
Col. D, 
table 2 


Oct. 
Col. E. 
table 2 


Nov. 
Col. D, 
table 2 


Nov. 
Col. E, 
table 2 


Dec. 
Col. D, 
table 2 


Dec. 
Col. E, 




Mar. 


Apr. 


May 


June 


July 


Aug. 


Sep. 


table 2 


For low/ 

For mean/. 

For high/_ 


10.9 
10.6 
11.7 


13.3 
12.8 
14.1 


13.7 
13.1 
13.9 


11.1 
10.6 
11.3 


14.3 
13.8 
14.4 


16.2 
15.5 
16.2 


23.5 
22.1 
23.8 


18.1 
17.1 
17.7 


17.8 
16.8 
17.5 


19.3 
18.3 
19.3 


13.0 
12.4 
13.0 


17.8 
17.0 
18.0 


16.1 
15.4 
16.3 


12.8 
12.3 
13.0 


12.4 
11.9 
12.6 


14.5 
13.8 
14.6 


12.2 

n.6 

12.3 



rnis deviation of prediction given in table 1 of 
McNish and Lincoln [1954], in which the deviations 
range from 7.0 to 19.2. Although these are of the 
same order of magnitude, it must be remembered 
that the McNish-Lincoln formulas were obtained 
using data from 10 sohir cycles, whereas the starulard 
errors of this ])a])er were cakmlated using less tlian 
2 cycles. Essentially the sunspot numbers used in 
the McNish-Lincoln method form a series of numbers 
spaced one solar cycle or about 11 years apart, 
whereas in (his method the sunspot numbers are 
spaced 1 Acnr apail . 

5. Results and Conclusions 

Since tliis paper was begun, values of maximum 
hourly median/o/'2 at Washington have been received 
for January tlirough May 1961. These are listed in 
table 5, tog(4her with the values of 12-month running 
average Zurich sunspot number estimated by the 
method of this paper, and those predicted by 
the McNish-Lincoln method. It must be remem- 
bered, of course, that the latter predictions were 
made, for lh(^ correspondhig months, 6 months 
earlier than tlu^ predictions using the critical fre- 
quencies at Washington. 

At the time of final version of this paper (June 
5, 1961) the provisional 12-month running average 
Ztirich sunspot number, centered on December 1, 
1960, had dropped to 84. :i. At the average rate of 
decline for the past 4 months, the number centered 
on the middle of January 1961 would be 76.7, 
compared witli tlie prtMlicted values 75.0 and 97.4 
of table 5. 



Table 5. Values of f at Washington, estimated values of R, 
predicted values of R by McNish-Lincoln method, Jan- 
Maij 1961 



1961 


/(Mc/s) 


R (Est.) 


R (Pred., McNish-Lincoln) 


Jan 

Feb 


9.2 
8.6 
8.6 
7.4 
6.8 


a 75.0 
59.4 
67.1 
59.5 
59.5 


97.4 
93.8 


Mar 

Apr -_ 


90.2 
88.5 


May 


86.2 







>* From table 2, col. 1); all others from table 2, col. P]. 

A quick method of this type sliould prove espe- 
cially valuable during the first 2 years of the rising 
j)art of a solar cycle, to help give an early indication 
of the level and timing of the coming maximum. 

It would of course be desirable to augment the 
sample size for estimating the regression equations of 
table 2 as more 3^ears of data become available, in 
paj'ticular the 3^ears 1961 and 1962, thus completing 
a 22-year period. Also it should be noted that a 
similar thing could be done for any ionosphere 
sounding station having a long run of well-taken 
data. Thus delays in waiting for information from 
other parts of the world could be avoided. 

An indication of the degree of accuracy at certain 
otiier locations may be obtained from table 6. Esti- 
mated values of R for January 1961 from the regres- 
sion equations given in this table, are, for tlic respec- 
tive stations, 70.6, 75.7, 81.4, and 84.5. 

It is suggested that adequate I'esults might be ob- 
tained from two locations, one in eacli hemisphere, 
at about the geomagnetic latitude of Washington, 
each station to cover its own winter months. For the 



Table G. Correlation coefficients, regression equations, standard errors for four stations, 1945 to 1960, Jan., J ulij 



Station 


Geomag. 
lat. 


Corr. coelT. 


Regression equations 


h Standard errors 




Jan. 


July 


Jan. 


July 


Jan. 


July 


Fairbanks 


64° N 

50° N 

48° S 

1° S 


0.98 
.99 
.94 
.97 


0.91 
.92 

a. 96 
.96 


/? = 96.04-21.2(/-8.3) 
/e=96.0+20.3(/-10.2) 
/^ = 96.0+48.7(/-7.7) 
/?=96.0+28.7(/-11.2) 


/?=98.5+64.4(/-5.7) 
/2=98.5+60.7(/-6.7) 
a/?=102.4+22.8(/-8.9) 
/?=98.5+35.7(/-8.9) 


13.9 
9.5 
22.2 
15.8 


27.7 


Washington 


24.3 


Christchurch 


'» 18.9 


Huancayo __ .- 


17.3 







a 15 years; data missing for July 1955. 

b Standard error of an estimated R for mean/. 
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equinoxes, possibly the station with the longer time 
series would be used. 

Finally, the presentation of two sets of equations 
in table 2 (columns D and E) is intended to be illus- 
trative of two ways of obtaining the result. One 
should not necessarily switch from one type of equa- 
tion to another, because in small samples such as 
these the reversal in size of correlation coefficients or 
other criteria could happen frequently and not be 
significant. 
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